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I. INTRODUCTION 



Black holes have been investigated in great depth and detail for more than forty 
years. However, almost all previous studies have focused on isolated black holes. On the 
other hand, one cannot rule out the important and more realistic situation in which black 
holes are actually embedded in the background of universe. Therefore black holes in non-flat 
backgrounds form an important topic. 

As early as in 1933, McVittie [1] found his celebrated metric for a mass-particle in the ex- 
panding universe. This metric gives us an concrete example for a black hole in the non-flat 
background. It is just the Schwarzschild black hole which is embedded in the Friedman- 
Robertson- Walker universe although there was no the notion of black hole at that time. In 
1993, the multi- black hole solution in the background of de Sitter universe was discovered 
by Kastor and Traschen [2]. The Kastor-Traschen solution describes the dynamical system 
of arbitrary number of extreme Reissner-Nordstrom black holes in the background of de 
Sitter universe. In 1999, Shiromizu and Gen extended it into the spinning version [3]. In 
2000, Nayak etc. [4,5] studied the solutions for the Schwarzschild and Kerr black holes in 
the background of the Einstein universe. 

In this letter, we extend the McVittie's solution into charged black holes. We first 
deduce the metric for a Reissner-Nordstrom black hole in the expanding universe; several 
special cases of our solution are exactly the same as some solutions discovered previously. In 
the previous work [6] we have applied the asymptotic conditions to derive the Schwarzschild 
metric in the expanding universe, which is exactly the same as that derived by McVittie 
by solving the full Einstein equations. That demonstrates the power of this simple and 
straight-forward approach. In this work we follow the same procedure to derive the metric 
for the Reissner-Nordstrom Black Holes in Friedman-Robertson- Walker Universe. We then 
study the influences of the evolution of the universe on the size of the black hole. Finally in 
order to study the motion of the planet, we rewrite the metric from the cosmic coordinates 
system to the Schwarzschild coordinates system. 
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II. DERIVATION OF THE METRIC 



The metric of Reissner- Nordstrom black hole in the Schwarzschild coordinates system 
is given by 



' ^ ' ' ' ^ ' ^ ' df^^f^ (de^ + sin^ edcf)^) , (1) 



where M and Q are the mass and charge of the black hole, respectively. 

For our purpose we rewrite the metric Eq.(l)in the isotropic spherical coordinates. We 
assume that = v and = x. So make variables transformation 



t = 2v, 5 = 2/, 2f = x(l + — ] - 



x J x'^ 



(2) 



then we can rewrite Eq.(l) as follows 



dl' 



(l-Ml + Q 



2\2 



-^dv + 



X J X'^ 



(dx^ + x^de^ + x" sin^ ed(l>'') . (3) 



_(l + f) - 

As is known, the metric for the FRW (Friedman-Robertson- Walker) universe is given by 



dF = -dv^ + 



"'^''^ , (dx^ + x^dd^ + x^ sin^ 



(4) 



(l + A;xV4) 

where a(v) is the scale factor of the universe and k gives the curvature of space-time as a 
whole. 

Taking account of equations Eq.(3) and Eq.(4), we set the metric for a Reissner- 
Nordstrom black hole embedded in the FRW universe as follows 



dl'^ = -A^ {v, x) dv"^ + B"^ {v, x) {dx^ + x^dO"^ + x^ sin 



2 /)^J,2 



(5) 



Then from equation Gqi = one obtains 



A{y,x)=f{y) 



B 

2B' 



(6) 
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where " • " denotes the derivative with respect to v. 

Compare the gu terms in Eq.(3) and Eq.(5), the possible form for the function B{v,x) 

is 

s{v) 



B (v, x) 



w (v, x) + 



X 



x^ 



(7) 



We note that the mass and charge of the black hole is concentrated in the singularity. In 
other words, there is no space distribution for mass and charge. Thus q and s which are 
related to the mass and charge, respectively, are only the functions of time v. 
Inserting Eq.(7) into Eq.(6), we obtain 

A.(v x) — V ^ ^' w'x w-'x-' 2w''x-' (^\ 

\ WX ) W^X'^ 

In the case oi v = const and the asymptotically fiat conditions, A{v = const, x) should 
be reduced to the yf—gm term in Eq.(3). Thus comparing Eq.(8) with the \/—gm term in 
Eq.(3), we infer the following identities should always hold 

wf 



1, 



w 



f 

(wq + wq) = 0, 



w'^x'^ \wxj 
sf s 



IllP'X^ iiP'x'^ ' 



(9) 



namely 



wf = w, qf = -q, sf = -2s, (10) 



Prom Eqs.(lO) we obtain 

^l + kxyi b h 62 

where M and Q are two integration constant which are related to the mass and charge of the 
black hole; h{v) is an arbitrary function which is related to the scale factor of the universe; 
the form of uj is obtained by inspecting the McVittie's solution. 



(dx^ + x^de"^ + sin^ Odcj)'^) , 



X 



(12) 



where a = a{v). 

Substituting Eqs.(7, 8, 11) in Eq.(5) , we obtain our final metric for tlie Reissner- 
Nordstrom black hole in the background of FRW universe 



di' 



1 - 



M2 



l + kx^/A)] +^(1 + ^x74) 



dv' 



(l + SVl + ^^V4) -£.(1 + ^x74) 

i+^;rT^)^^(i+..V4) 



{l + kx'^/Ay 
(dx^ + x^de"^ + x^ sin^ ed(i)'^) , 



(13) 



where we have made a variable replacement h{yY — > 0(1;). Eq.(13) is derived from Eq.(7) 
and when a{v) — constant, k — Eq.(13) restores Eq.(3). So Eq.(3) satisfies Eq.(7). 

For the Reissner-Nordstrom-de Sitter metric, we have a — e^*, A; = 0, so Eq.(13) becomes 
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dV 



I 

~^ — 5" ~r ~i — 7 



\ ax J 



;dv + a 



ax J 



a^x"^ 



{dx^ + x^de"^ + x^ sin^ 



(14) 



We will show in the next section that Eq.(14) can be reduced to the familiar form in 
Schwarzschild coordinates. When H = 0, Eq.(14) restores the Reissner-Nordstrom met- 
ric as given by Eq.(3). When Q = 0, Eq.(13) becomes 



dl^^ 



1 - 



M 



1 + fcxV4 



rdv'^ 



(l + M^l + A;xV4)^ (l + A;a;V4)^ 

(dx'^ + x'^de'^ + x^ sin^ Odcjp) , 



1 + — x/l + A;xV4 
ax ^ 



(15) 



It is just the McVittie solution. Another special case of our solution is for the extreme 
Reissner-Nordstrom black hole, M = Q, in the de Sitter universe. In this case, Eq.(14) is 
reduced to a special case of the Kastor and Traschen solution [2] for a single black hole. 



(i + m 

\ ax J 



ax 



(16) 



In Eq.(16), a = e^*. It describes one extreme Reissner-Nordstrom black hole in the de Sitter 
universe. 

III. FURTHER DISCUSSION ON THE METRIC 



In this section we will verify that Eq.(13) satisfies Einstein-Maxwell equations. The 
Einstein-Maxwell equations may be written as 

Gnu — Sn {Thj, + E^j,) , 



F^r = 0, 



A —A 



(17) 



where TJ^i^ and E^j, are the energy momentum for the perfect fiuid and electromagnetic fields, 
respectively, which are defined by 



E = — (f F"--o F aF°'^ \ 



(18) 



where p and p are the energy density and pressure. 11^ is the 4- velocity of the particles. F^^/ 
and are the tensor and the potential for electromagnetic fields. 

Input the components of the metric, Eq.(13), to the Maple software package, we obtain 
the Einstein tensor G^jy and the energy momentum tensor T^^ and E^^ for the perfect fluid 
and the electromagnetic fields 



^0 — Pi 
I'kEq — 1'kE\ — 

2ttEo — 2ttE\ 



-'l — -'-2 — -'S ~ Pi 



g2 (1 + kx^/A) 



(l + SVl + ^^V4) -^(1 + ^x74) 
Q2 (1 + kx^/A) 



4' 



x'^a* 



(1 



M 
ax 



l4- 



l + kx^/A) -^{l + kxyA) 



(19) 
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Substituting the above components of electromagnetic tensor in the second equation of 
Eqs.(18), we obtain the non- vanishing components of electromagnetic tensor F^^, 



1 - (1 + kx^l^) + ^ (1 + kx^l^) 

^ (20) 



Then substituting Eq.(20) in the second equation of Eqs.(17), we obtain the non- vanishing 
components of the potential 

Ao = I F'^'googudx. (21) 

It is a straightforward work to verify that Eq.(20) also satisfies the last equation in Eqs.(17). 
Since = always holds, thus from Einstein equations we have = T!^'^ + Et",^ . On the 
other hand, we have the relation 

AnEtl'^ = F>^^F^.^ + F.lt'^F^ - ]^g^^F^PF^p 

~ ^ a;iy ' 2^ ^ P'^'^ P^''^ '^'^P' 

= -Fi^^Jo, = 0. (22) 

So both T^y and E^^ satisfy Bianchi identity. We therefore conclude that Eq.(13) is an exact 
solution of the Einstein- Maxwell equations. 

To show how the two parameters M and Q are related to the mass and charge of the 
black hole, we assume the evolution of the universe is much slower and approximately adopt 
the mass formula for the stationary space-time [8] (To our knowledge, we can only define 
the mass for the stationary and asymptotically flat space-time) . We find the mass Mq and 
charge Qo of the black hole are given by 

Mo = / ea.^VX" = — , 
STT Js a 



47r Js 



^abcdF'^'^ — —■ (23) 

s a 

Thus for the observer in the infinity, the black hole's mass and charge will decrease with the 
expansion of the universe and increase with the contraction of the universe. We will return 



to this point in the following discussion again. 

Let's now consider the two typical surfaces of the black hole in the cosmic coordinates 
system. We obtain the radius of the time-like limit surface (TLS) [8] 



rTLS = J—-^: 24 



and the time derivative of the radius of the event horizon [8] 



1 , 

J- — 9 5 — ~r 9 5 



Teh = ±- 



2 



M Y _ _Q2_ 



2 ' 



(25) 



where two signs " + " and " — " correspond to the expanding and the contracting universe, 
respectively. Since the event horizon is always in the inner of TLS, we have teh < ^ for 
expanding universe and teh > for contracting universe. 

Eqs. (24, 25) tell us the typical scales of the black hole are closely related to the evolution 
of the universe. They shrink with the expansion of the universe and expand with the 
contracting of the universe. For asymptotically flat background, a — 1 and fEH — 0, these 
two kind of surfaces coincide 



VTLS - TEH = - Q'- (26) 

Comparing Eq.(26) with Eq.(24), we find that the black hole in the expanding universe has 
the mass Mq = M/a and charge Qo/a. They are both dependent on the scale of the universe 
and not a constant. 

Prom Eq.(23), one obtain the lifetime of a black hole 

r^^. (27) 
It is approximately the age of the universe. 



IV. EQUATION OF MOTION OF A PLANET 



In this section, we turn to the discussion of the problem of the orbit for a planet 
in the expanding universe. Current measurements of the microwave background radiation 
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show that our universe is highly hkely flat in space [7] . So in the next we consider the local 
dynamics in the space-flat space case, ie., k — 0. Eq.(13) is then reduced to 



\ ax J a'^x-' 



1 2 



dv + a 



ax J a'^x^ 



(dx^ + x^dO"^ + x^ sin^ Odct)'^) 



(28) 



In order to study the motion of a planet, we should rewrite Eq.(28) in the Schwarzschild 
or solar coordinates system. Similar to Eq.(2), make variables transformation as follows 



M\2 

T^2v, s^2l, 2r^ax(l + — ] - 



ax J a'^x'^ 



(29) 



Then Eq.(28) becomes 



ds' 



( 2M g2 . o\ , o , 2M Q2 2 
- 1 + ^ - MT" + 1 + ^ "^^ 



( 2M Q 
-2rH 1 + 



2\ -1/2 



dTdr + {de"^ + sin^ Odcj)^) , 



(30) 



where H = is the Hubble parameter. The coordinate system of {T,r,9,(l)) is not 
orthogonal. We can eliminate the coefficient of dTdr by introducing a new time coordinate, 

t. The form of Eq.(30) suggests we set 



. X . X ( 2M Q 
dt^F (T, r)dT + F (T, r)rH 1 + 



[l- — + ^-H\'\ dr, (31) 



where F{T,r) is a perfect differential factor and it always exists. Then Eq.(29) can be 
written as 



ds^^ 



2M 



F^dt^ + I 1 



2M 



+ ^-H'r'\ dr' + r'dn\ (32) 



where H and F are both the functions of variables t and r. If p = p = 0, we have 
H — 0. Eq.(30) and Eq.(32) both turn into the static Reissner- Nordstrom solution. If 
2M/r — )• and jr^ — > 0, Eq.(32) represents the solution for the Friedman- Robertson- 
Walker universe in Schwarzschild or solar coordinate system. For the de Sitter universe, H is 
a constant. Eq.(31) tells us we may choose F(i, r) — 1. Then Eq.(32) is just the well-known 
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Reissner-Nordstrom-de Sitter space-time. It is a static space-time. Thus the geodesic is 
time-independent. In other words, the orbit of the planet does not vary in the background 
of de Sitter universe. 

In general, H is the function of t and r . So Eq.(32) is a non-stationary space-time and 
the geodesic in this space-time is generally time-dependent. Thus we conclude that universe 
expansion would influence the orbit of planet. However, the effect term H^r^ is very small. 
Since 

. _ H'r' _ P 9 
~ 2M/r M/(47rrV3) ~p ^ ^ 

where p is the universe energy density and p is the energy density within the region r where 
the energy M is distributed. The ratio of p to p varies from 4 x 10~^^ for Mercury-Sun 
system (In detail, r represents the distance from the Sun to Mercury and M is the mass of 
the Sun. The mass of the planet is much smaller than the Sun.), 1.8 x 10~^^ for Neptune-Sun 
system to 10~^ for Galaxy. 

Keep this in mind or regard H as a constant, we obtain the motion equation of a planet 

+ " = ^ + 3M.^ - 2Q\^ - (34) 

where u = 1/r and the prime denotes differentiation with respect to 0. L is the angular 
momentum of the planet. The term M/ Li^ is the most significant one. The two terms 3Mm^ 
and H'^/{LF'u^) come from the GR (general relativity) effect and CE (cosmic expansion) 
effect, respectively. Q'^ujl? and —IQ'^u^ are related to the charge of the source. Eq.(34) 
tells us the orbit of the test body will be influenced by the expansion of the universe. Since 

^^^^ ' = 2£ = 8 X 10-^^ for Mercury, 3.6 x lO^^^ for Neptune, (35) 
so the influence is related to the ratio of p to p. It is very small and neghgible. 



V. CONCLUSION AND DISCUSSION 



In conclusion, we have presented the metric for the Reissner-Nordstrom black hole 
in the background of FRW universe. It extends McVittie's solution, Reissner-Nordstrom- 
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de Sitter solution and the special case for a single black hole of the Kastor and Traschen 
solution. Assume the evolution of the universe is much slower and adopt the formulas for 
computing the mass and charge of stationary space-time, we find that both the mass and 
charge of the black hole decrease with the expansion of the universe and increase with the 
contraction of the universe. We also find that the two typical scales, the time-like surface 
and the event horizon, of the black hole both shrink with the expansion of the universe and 
expand with the contraction of the universe. This is due to the fact that the mass and 
charge of the black hole are both varying with the evolution of the universe. 

To obtain the equation of motion of a planet, we rewrote the metric from the cosmic 
coordinate system to the Schwarzschild or solar coordinate system and deduced the geodesic 
equation. The equation shows that the orbit of the planet in the Reissner-Nordstrom field 
embedded in the FRW universe will be influenced by the evolution of the universe. The 
magnitude of influence depends on the ratio e between the energy density of the system 
and the energy density of the universe. Since the ratio e is extremely small, varying from 
4 X 10~^^ for Mercury-Sun system, 1.8 x 10~^^ for Neptune-Sun system to 10~^ for Galaxy, 
the influence of the expansion of the universe is very small and negligible [9] . 
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